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ABSTRACT 

We present an analysis of approximately 200 hours of observations of the pulsars 
J1634—5107, J1717—4054 and J1853+0505, taken over the course of 14.7 yr. We show 
that all of these objects exhibit long term nulls and radio-emitting phases (i.e. minutes to 
many hours), as well as considerable nulling fractions (NFs) in the range ^ 67 % — 90 %. 
PSR J1717—4054 is also found to exhibit short timescale nulls (1 — 40 P) and burst 
phases (< 200 P) during its radio-emitting phases. This behaviour acts to modulate the 
NF, and therefore the detection rate of the source, over timescales of minutes. Furthermore, 
PSR J1853-I-0505 is shown to exhibit a weak emission state, in addition to its strong and 
null states, after sufficient pulse integration. This further indicates that nulls may often only 
represent transitions to weaker emission states which are below the sensitivity thresholds of 
particular observing systems. In addition, we detected a peak-to-peak variation of 33 ± 1 % 
in the spin-down rate of PSR J1717—4054, over timescales of hundreds of days. However, no 
long-term correlation with emission variation was found. 

Key words: pulsars: general — pulsars: individual (J1634—5107, J1717—4054, 

J1853H-0505). 


1 INTRODUCTION 


Pulsars are classically considered to be rapidly rotating neutron 
stars, which radiate beamed radio emission from their poles in a 
steady and predictable fashion. In this model, an observer on Earth 
receives a pulse every rotation from a given pulsar, as its beam of 
emission sweeps across the observer’s line-of-sight (LOS). In re¬ 
ality, however, pulsars have been shown to exhibit variability in 
their emission over every timescale which they can be observed 
(i.e. from nanosecond bursts to multi-decadal vari ations; see, e.g., 
iHankins et 'Zll2003l:[K^l2013l:lLvne ^l2013h . 

Of particular interest is the phenomenon of pulse nulling, 
where the pulsed em ission from a pulsar appears to completely 
cease teackeJ [I?^ . over timescales of pulse periods ( P) to 
many years (see, e.g.. lWang et al.ll2007l : lLorimer et alj|2012h . This 
can be caused by a number of scenarios: 1) a pulsar could un¬ 
dergo complete ce ssation (a.k.a. deep nulls; e.g. lKramer et alj2006l : 
lOaiiar et au 120121) or 2) transition to a weaker e mission mode, 
which is below a given sensitivity threshold (e.g. lEsamdin et al.l 
l2005LlYoung et al.l20l4) : 3) the radio beam could also move out of 
the LOS due to unfavourable emission geometry (e.g. Dvks et afl 
I 2 OO 5 L iMelikidze & Giill2006l: IZhang et Ixhmi: lTimokhinll2o7^ 
or 4) the acceleration zone might not be completely filled by 
electron-positron pairs, resulting in time-dependent variations in an 
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emission ‘carousel model’ (see, e.g., IPeshnande & Ran^l200ll : 

Ijanssen & van Leeuwenll2004:lRanidn & Wrighlll2007). 

While the latter geometrical and beam-filling theories can 
account for short nulls (i.e. < 10 P), they cannot account for 
the typical (~ 10^“^ P) to long-term (~ 
serve d in many objects (e. g. iKramer et al.ll200^ : iMcLaughlin et al.l 
I 2 OO 6 I : ICamilo et al.l l2Q12h . In this context, nulls are considered 
to be extreme manifestations of mode changes in p ulsar emis¬ 
sion, due to global m agnetospheric state changes dWang et'H] 
l2007l : lLvne et al.f2010l) . This is strongly supported by the observed 
correlations between emission variability and spin-down {v) rate 
changes in several sources (e.g. iLvne et al.l2010l : TLi et al.l2012h . as 
well as the simul taneous X-ray—radio mode switches observed in 
PSR B0943-I-10 dHermsen et al.ll2013h . However, the triggers re¬ 


sponsible for such magnetos pheric reconfigurations have still not 
been firmly ascer t ained (e.g. Icordes & Shannonll2008l : lUvne et ^ 

l2010l:lj one 31201 ll:lRosenetal.ll201ll). 

Significant emphasis has been placed on discovering and 
monitoring pulsars with extreme nulling and/or moding proper¬ 
ties, which may give valuable insight to determine the root cause 
of pulsar emission variability. Such objects may undergo mode 
changes and/or nulls with long timescales (e. g.lKram er et alj2006l: 


I Keith et alj2013l : [Knispel et alj2013l : lBrook et alJ gOlj ). They may 

also exhibit large nullin g fractions (NEs: e.g. lBurke-Snolaor et al.l 
l201ll : l^aiiar et al.]|2012h . However, very few sources are known to 
display such features (see lBurke-Snolaor et alJl^l ll for a detailed 
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discussion), which has somewhat stifled the amount of studies on 
these enigmatic objects. Below, we shift our focus to three pulsars 
which have been shown to exhibit extreme nulling behaviour. 


The first pulsar we discuss, J1634—5107, is a 507 ms source 
which was discovered in th e Parkes Multibeam Pulsar Survey 
(PMPS: iLorimer et alJliooeh . Nulling activity in the source was 
first identified in further analysis of a number of PMPS discov¬ 
eries, where it was shown that the pulsar exhibits a strong radio- 
emitting (i.e. ‘ON’) state and a pure n ull (i.e. ‘OFF’) stat e with 
a periodicity of approximately 10 days jO’Brien et alJl2006h . This 
analysis was later elaborated on bv lO’BrienT i 2010h . who analysed 
many hours of observations — recorded over hundreds of typically 
short (< 30 min) sessions — to infer a NF ~ 86 % for the source. 
Further investigation of the emission phases also showed that the 
object does not fluctuate between null and emitting states on short 
timescales (i.e. fon < 1 d and tos < 9 d). 


PSR J1717—4054, the second pulsar, is a 888 ms source 
which was discov ered in a Southern Ga lactic plane survey with the 
Pa rkes telescope J johnston et al.|[l992h . Follow-up timing studies 
bv ijohnston et al.l ( Il995h showed that the source undergoes pulse 
nulling for ~ 75 % of the time. Consequently, the number of detec¬ 
tions were too few to e nable a measuremen t. This source was also 
another subject of the lO’Brien et al.l (l2006h PMPS study, where a 
slightly higher NF ~ 80 % was inferred. In this study, it was noted 
that the source transitions to and from null phases over timescales 
less than a second, indicating an abru pt cessation mechanism. Fol¬ 
lowing this work, IWang et ^ ( l2007h presented the analysis of a 
single 2-h observation of the source. They detected the object for 
the first 210 s of the observation, after which the object transi¬ 
tioned to, and remained in, its null state. As a result, they inferred 
a NF > 95 % and transition timescale — the expected length of 
time between consecuti ve mode t ransiti ons — of > 2 h. This re¬ 
sult was later refined by lO’BrienI ( l2010ll . who inferred an average 
NF ~ 74 % from the analysis of 252 observations. They also re¬ 
ported the presence of short timescale nulls during active emission 
phases. 


The last pulsar, J1853-I-0505, is a relati vely undocumented 
905 ms source. It was discovered in the PMPS dHobbs et al.ll2004h 
and, since, has only been documented fo r its pulse broad ening 
characteristics due to interstellar scattering dBhat et al.ll2004h . The 
nulling activity in this source was only identified through further 
analysis of the PMPS discoveries and has remained unpublished 
until now. 


Since their initial discovery, PSRs J1634—5107, J1717—4054 
and J18534-0505 have been routinely observed under several ded¬ 
icated observing programmes, predominantly with the Parkes tele¬ 
scope. This has facilitated more comprehensive analysis of their 
emission variability and timing pro perties, whi c h is p resented in 
this work. We note, however, that iKerr et ai] d2014ll published 
results from an independent analysis of the same data set for 
PSR J1717—4054, during the preparation of this manuscript. 

In the following section, we review the observing programmes 
used to gather our data. We then discuss the emission, timing and 
polarimetric properties of the objects, where possible, in sections 
mil an dm respectively. Lastly, we discuss the implicati ons of our 
results in section[6]and compare them with the results of iKerr et al.l 
(I 2 OI 4 I) for PSRJ1717-4054. 


2 OBSERVATIONS 

Observations of PSRs J1634—5107 and J1717—4054 were ob¬ 
tained through an intermittent source monitoring programme (IMP) 
which was conducted with the Parkes 64-m radio telescope, over 
the period of 1999 August 21 to 2014 May 12. The majority of 
these data were obtained with the H-OH receiver and central beam 
of the Multibeam receiver. However, a number of observations were 
also obtained using the 10/50cm receiver (refer to TableQ] see also 
iManchester et al.ll2013l for detailed specifications on the Parkes re¬ 
ceivers). 

An analogue filterbank (PAFB) system was primarily used to 
record the observations for these pulsars up until 2008 September, 
after which digital filterbank (PDFB) systems were used to obtain 
the majority of the remaining data (see Table [T}. Note that reg¬ 
ular observations of the pulsars with the PDFBs commenced from 
2007 December. A total of 14 observations were also recorded with 
the ATNF Parkes Swinburne Recorder, CASPER Parkes Swinburne 
Recorder and Wide Band Analogue Correlator backends for these 
pulsar^ 

The PAFB observations were one-bit digitised at 250 fj,s inter¬ 
vals and were later folded offline to produce both folded data, with 
typically 256 bins per period at 1 min sub-integration intervals, as 
well as single-pulse archives where possibl^B The PDFB observa¬ 
tions were folded online to typically provide 512 bins per period, 
using 8-bit digitisation, at 1 min sub-integration intervals. A po¬ 
larised calibration signal was also injected into the receiver probes, 
and observed, prior to a large number of these observations with 
the Multibeam receiver, in order to polarisation calibrate a subset 
of the data. 

Observations of PSR 118534-0505 were obtained with both 
the Parkes 64-m and Lovell 76-m radio telescopes, over the period 
of 2000 March 30 to 2014 May 12 under a joint IMP. The same 
receivers and observing setup as described above was used to ob¬ 
tain observations of this source at Parkes. The Lovell observations 
were obtained with a 20-cm dual-orthogonal, linear feed receiver, 
coupled with an analogue filterbank (LAFB) or digital filterbank 
(LDFB). The LAFB data were folded online to provide 400 bins 
per period at 1 min sub-integration intervals, using one-bit digitisa¬ 
tion. The LDFB observations were folded online to provide 1024 
bins per period at 1 min sub-integration intervals, using 8-bit digi¬ 
tisation (refer to Table[T}. 

In off-line processing, we de-dispersed and examined the data 
for radio-frequency interference (RFI). For the Parkes and LDFB 
data, we carried out manual RFI mitigation through use of PSRZAP 
and PA^ so as to reduce the number of frequency channels and sin¬ 
gle pulses/sub-integrations excised from the data. We also reduced 
the number of bins to 256 per period for each of these observa¬ 
tions, for the emission variability analysis. Whereas, we reduced 
the number of bins to 150 per period for the polarimetric analysis. 

For the LAFB data, we were restricted to only 32 MHz of 
bandwidth. Therefore, limited manual RFI excision was only possi- 


^ Refer to http://astronomv. swinbume.edu.au/pulsar/?topic=instmmentation 
and iManchester et al.l i2013l) for more details on Parkes backends. 

^ Single-pulse data was obtained for 331 observations of 
PSR J1634-5107, 205 observations of PSR J1717-4054 and 90 ob- 
servati ons of PSR J1853-4 0505, respectively. 

® See iHotan et al.l j2004l and http://psrchive.sourceforge.net/manuals for 
details on the PSRCHIVE software suite. 
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Table 1. The observation characteristics of the observing programmes conducted for each source. PSR denotes the pulsar observed, REC denotes the receiver/ 
backend reference, MJD refers to the modified Julian Date at the start of the observations and Tapan represents the total time-span of the observations. The 
total number of observations carried out is given by Af^ba, T^ba denotes the typical observation duration and (C) represents the average observation cadence. 
The centre sky frequency, observation bandwidth and total number of frequency channels are also denoted by u, Au and Afghan, respectively. 


PSR 

REC 

MJD 

Tspan (d) 

Nobs 

Tobs (min) 

(C) (week 

u (MHz) 

Av (MHz) 

'^chan 

J1634-5107 

Multi 

51411.4 

5378.2 

383 

7.6 

0.5 

1374 

288 

96 


H-OH 

52961.1 

281.2 

73 

10.3 

1.8 

1518 

576 

96 


10/50cm 

53000.2 

609.0 

35 

10.7 

0.4 

2934/685 

576/64 

192/256 

J1717-4054 

Multi 

53280.3 

3509 

335 

5.5 

0.7 

1374 

288 

96 


H-OH 

54112.0 

116 

21 

14.3 

1.3 

1518 

576 

1024 


10/50cm 

53029.8 

3181 

15 

5.0 

0.03 

3094/732 

1024/64 

1024 

J1853-I-0505 

Multi 

51633.8 

3212.8 

94 

8.6 

0.2 

1374 

288 

96 


H-OH 

53150.6 

90.9 

19 

14.5 

1.5 

1518 

576 

1024 


lOcm 

53190.5 

1.2 

6 

15.0 

35.0 

685 

64 

1024 


LAFB 

51843.4 

3527.7 

250 

16.5 

0.5 

1396 

32 

32 


LDFB 

54846.6 

829.7 

23 

11.8 

0.2 

1520 

384 

768 


ble for these data. The LAFB data were also converted into the SIG- 
PROC0 file format, to facilitate variability analysis of the source. 


3 EMISSION PROPERTIES 

3.1 Variability and Emission Timescales 

Here we describe the nulling activity of the three pulsars chosen for 
our study. For each of these objects, the inferred s cintillation band¬ 
width from the NE2001 model is less than 2 kHz dCordes & Lazid 
l2002h . This is consistent with our analysis of their dynamic spec¬ 
tra, which do not exhibit any prominent scintles given the frequency 
resolution afforded (i.e. > 500 kHz). Thus, the variability of these 
objects cannot be described by propagation effects in the interstel¬ 
lar medium (e.g. scintillation). Rather, the nulling behaviour de¬ 
scribed below is considered to be intrinsic to the sources. 

3.1.1 PSRJ1634-5107 

PSR J1634—5107 was detected in 64 out of the 491 observations 
analysed, which amounts to a total of ~7.0 h in the ON state from 
the ~67.7 h data set. No weak or underlying emission was observed 
during individual OFF observations, or in the ave rage OFF pro¬ 
file fo r each observing frequency (c.f. B0826—34; lEsamdin et al.l 
l2005h . which is consistent with previous findings. Moreover, we 
discovered that the detectability of the source is strongly corre¬ 
lated with integration length. That is, we found that sub-integration 
lengths of < 2 min do not provide sufficient sensitivity for re¬ 
liable characterisation of mode transitions due to apparent nulls. 
This can be explain ed by the low flux density of the sour ce (c.f. 
PSRB1931+24; e.g. lKramer et al.l2006l:l^ung et alj2Qlil and/or 
short timescale nulling not probed by our data. 

Therefore, we only comment on the properties of the object 
on an observation to observation basis, where the mean length is 
~8 min. We subsequently estimate a NF = 90 ± 5 % for the 
source — from the fractional time it was obse rved in t h e OFF 
phas^B ~ which is consistent with the result of lO’Bri^ ( I 2 OIOI 1 . 
Note that we quote the uncertainty on the NF assuming Poisso- 
nian sampling statistics. That is, we assume that each measurement 

http://sigproc.sourceforge.net/ 

® Under the assumption that pulse nulling is a stochastic process, a NF can 
be accurately characterised from the fractional time observed in an OFF 
phase over a sufficient number of observations. 



Figure 1. Average emission profiles for PSR J1634—5107 from the longest 
observations, separated by emission mode and normalised by their peak 
intensities. Bottom: OFF profile from an 4.5 h observation taken on 
2002 September 13 at 04:52:05 UTC. Top: ON profile from a 33 min ob¬ 
servation recorded on 2002 January 19 at 19:12:26 UTC, offset vertically 
for clarity. Note that the flux density of the OFF emission is < 1.7 ± 0.5 % 
that of the ON emission for these data. 

of the source in a particular emission state is independent. Thus, 
we estimate the fractional error on the NF as ~ V-^obs/A^obs- We 
stress, however, that this theoretical estimate will not be as accu¬ 
rate as a constraint obtained from documenting a large number of 
emission phases continuously. 

Further to the above, we did not observe any mode transitions 
during individual observations of the source. As such, the maxi¬ 
mum observation durations of 4 h 30 min for the OFF mode and 
33 min for the ON mode (see Fig. |T]( do not accurately constrain 
the source’s emission timescales. Instead, we place constraints on 
the emission timescales through analysis of the longest and highest 
cadence observing runs. 

In our data set, there were two particularly long (i.e. 10 h) 
observing runs (between 2002 November 1 and 2002 Novem¬ 
ber 3) where the pulsar was undetected throughout each half-hour 
to hourly observation. By comparison, there was one particularly 
long (i.e. 6.5 h) observing run on 2004 May 25 where the source 
was detected throughout each half-hour to hourly observation. Fur¬ 
ther to this, we place tentative limits on the maximum emission 
timescales from a series of observations, between 2003 August 12 
and 2003 August 16, where the pulsar was seen to switch between 
modes on a daily basis. This indicates that the source undergoes 
relatively long emission phases, which are at least several hours up 
to a day in length. 
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3.1.2 PSRJ1717-4054 

PSR J1717—4054 was detected in 119 out of the 371 observations 
analysed, and was observed to transition between its ON and OFF 
modes on 38 occasions. Analysis of the sub-integration data for 
each observation showed that the source was ON for a total of 8.6 h 
of the 36.7 h data set. This equates to a NF = 77 ± 5 % — us¬ 
ing the error estimation method pres ented in ij 13.1.11 — w hich 
is consistent with th e findings of both Ijohnston et al.l ( Il995h and 
lO’Brien et al.l ( l2006h . We also infer the average transition timescale 
to be ~56 min from the total observation length and number of fully 
resolved mode transitions (Tobs/Aswitch). 

The null lengths in the sub-integration data were observed to 
span from 0.5 to > 64 min. Whereas, the ON phases were observed 
to last between 0.5 and 12.5 min in individual observations. Given 
the poor constraints on the upper limits of each emission timescale, 
we also analysed successive contiguous observations which were 
not separated by more than 10 min. From this analysis, we infer a 
maximum ON timescale of 16.2 min and the same maximum OFF 
timescale of > 64 min. Note that this analysis does not completely 
constrain the possible range of ON timescales due to inadequate 
observation cadence. Therefore, it is possible the maximum ON 
timescale could be longer than that inferred from this work. The 
variability of the source is demonstrated in an hour-long observa¬ 
tion shown in Fig.|2 

Upon inspection of the available single-pulse data, we were 
also able to confirm that the source transitions between emission 
modes over timescales of < IP. In addition, we found that the 
source exhibits a bimodal distribution of nulls. That is, the source 
preferentially undergoes nulls over timescales of ~ 1 — 40 P or 
> 340 P. We also note that emission bursts from the pulsar only 
last ~ 1 — 200 P. Therefore, it is clear that the majority of the short 
timescale nulls will not be resolved in the sub-integration data of 
the ON phases where fsub > 60 s. Thus, the interpretation of the 
ON emission timescales is dependent on the observation s etup. 

_ The above results are consistent with the results of lO’BrienI 

(l2010h. but are m arkedly different from the NF estimate of 
IWane et ^ ( l2007h taken from a single observation. We attribute 
this to variation in the OFF and ON timescales, where the wide 
range of activity timescales observed in our data set can lead to 
variation in NF estimates between observations. 


3.1.3 PSRJ1853+0505 

PSR J1853-1-0505 was detected in 93 out of the 392 observa¬ 
tions analysed, and was observed to transition between its ON and 
OFF states on 29 occasions in individual observing runs (see, e.g.. 
Fig. [3}. Analysis of the sub-integration data showed that the pul¬ 
sar was ON for a total of 17.8 h of the 92.8 h combined data set. 
For the LAFB data set we further integrated observations to sub¬ 
integration times of > 6 min, to counter the low bandwidth and 
subsequent reduction in sensitivity. Whereas, we were still able to 
probe typical timescales of > 1 min for the Parkes and LDFB data 
given the greater sensitivity. 

In addition to the strong detections, we noted that the source 
appears extremely weak in 12 of the observations (refer to Fig. 0. 
This weak emission is not detected over timescales of less than sev¬ 
eral minutes. Rather, it can only be observed through integration 
of sufficiently long observations (i.e. > 10 min), which is similar 
to th e underlying emission observed in a number of sources (see , 
e.g.. lEsa mdin et all 1 20051 : IWang et ^ l2007l : lYoung et^ 120141 : 
ISobev et alj2014h. 




Figure 2. An hour-long observation of PSR 11717—4054 recorded on 
2012 October 10 at a centre frequency of 3094 MHz. Top panel: Sub¬ 
integration versus phase intensity map showing the pulsar' transition be¬ 
tween the ON and OFF modes and back again. Bottom panel: Average 
profiles for the corresponding OFF (bottom) and ON (top) sub-integrations 
displayed in the top panel, normalised by their peak intensities and offset 
vertically for clarity. Note that no emission is detected through integration 
of 48 min of data. The null confusion limit is also < 0.8 ± 0.2 % of the flux 
density of the ON emission for these data. 


Remarkably, weak-mode transitions were also observed in two 
observations, spaced over a year apart. In the first observation, the 
pulsar exhibited weak-emission for ~ 16 min, then abruptly tran¬ 
sitioned to its strong mode (i.e. between ~ 1 min sub-integrations) 
for the remaining ~ 19 min of the observation. In the other ob¬ 
servation, the source emitted in its strong mode for < 1 min, then 
abruptly switched to its weak mode for the remaining ~ 9 min. 
These attributes indicate that the weak emission mode is relatively 
stable, which occurs between strong and OFF phases of emission. 
The above also suggests a common driving mechanism for the sep¬ 
arate emission modes, due to the consistency of their transition 
timescales. 

Although the LAFB data was subject to more RFI, compared 
with data obtained with the other backends, the longer observing 
runs recorded provide greater constraints on the cumulative ON and 
OFF timescales of the source. In these data, we did not resolve any 
ON modes of less than 30 min in duration, nor did we observe any 
greater than 2 h in length. Similarly, the minimum and maximum 
OFF timescales were found to be 40 min and > 7 h 20 min, re¬ 
spectively. However, we emphasise caution in assuming the upper 
limit for the maximum OFF length, given that the LAFB data were 
limited in being able to probe weak-mode detections. 
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Figure 3. A 34.4 min long observation of PSR J1853+0505 recorded 
on 2001 June 17 at a centre frequency of 1374 MHz. Top panel: Sub¬ 
integration versus phase intensity map showing the pulsar transition be¬ 
tween the OFF and ON modes. Note that the blank horizontal lines repre¬ 
sent sub-integrations weighted to zero from the RFI excision process. Bot¬ 
tom panel: Average profiles for the corresponding OFF (bottom) and ON 
(top) sub-integrations displayed in the top panel, normalised by their peak 
intensities and offset vertically for clarity. Note that no emission is detected 
through integration of 9 rain of OFF data. The null confusion limit is also 
< 3 ± 1 % of the ON emission for these data. 


We calculated NF estimates for both the LAFB (NFlapb = 
86 ±6 %) and combined Parkes-LDFB = 67 ± 8 %) data 

sets separately, due to the different observation sensitivities. Con¬ 
sidering the prevalence of weak emission in the combined Parkes 
and LDFB data set, i.e. ^11 % of the total time, we suspect that 
the NF estimate for the LAFB data is most likely overestimated. 
As such, we assume the NF from the combined Parkes-LDFB data 
provides a better representation of the variability of the source. We 
also infer an average timescale of ~3.2 h between contiguous ON 
modes from the total observation length and number of observed 
transitions. 


3.2 Flux Density Limits 

The detection of particularly weak emission in a number of ob- 

[ jects — through sufficient pulse integration ( see, e.g.. lEsamdin et all 
ZOO.'iL IWang et alJl2007 : lYoung et alJl 2014l) or migration to lower 
observing frequencies l Sobev et al.ll2014l) — raises speculation as 
to whether all nulls actually represent weak emission states or not. 
This is punctuated by the fact that every telescope has a specific 
flux density limit, for a given integration time, which may or may 
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Figure 4. Two 15 min observations of PSR J1853-I-0505, obtained at 
1518 MHz, demonstrating the different emission states of the source. The 
top plot shows a weak detection of the source, which was obtained an hour 
before an OFF detection and 2 hours before the strong detection, shown in 
the bottom plot. The integrated profile and frequency versus phase intensity 
plot is displayed for each observation, in the top and bottom panels of each 
plot, respectively. 


not allow the detection of such low intensity emission states. With 
the above in mind, we sought to place upper limits on the flux den¬ 
sities of the separate emission modes (S'ofi : Son) for the pulsars in 
this study, so as to characterise the possibility of null confusion in 
our data. 

To estimate the flux density of the sources, during their sep¬ 
arate emission modes, we formed time- and frequency-averaged 
profiles for corresponding observations obtained with the Parkes 
telescope. Observations were aligned using the timing solution for 
each source, prior to integration. We then used the modified ra¬ 
diometer equation to estimate the average flux densities attributed 
to each mode from the me an signal-to-noise ratios (SNRs; see, e.g., 
iLorimer & Krame^l2005t) . Here, we assume that the dominant un¬ 
certainties arise from gain and system temperature variations with 
respect to elevation for a given observing session. We thus con¬ 
servatively assume a 20 % variation in flux density between each 
observation. We note, however, that the actual uncertainites on our 
flux measurements are likely to be lower than this upper limit, espe¬ 
cially for the 10/50cm observations which were obtained in close 
succession. Table shows the result of this analysis for the three 
sources. 

Both PSRs J1634—5107 and J1853+0505 were not observed 
during their ON modes while the 10/50cm receiver was in use. 
Therefore, we cannot provide any constraints on the ON or OFF 
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Table 2. Average pulse properties for PSRs J1634—5107, J1717—4054 and J1853+0505. The pulsar J2000 designation, observation frequency, bandwidth, 
number of observations used for each fold and total integration length are denoted by PSR J, u, Au, A^on,off ^on,ofF’ respectively. The mean signal-to- 
noise ratio, average equivalent width and flux density are represented by SNR, Weq and S^n off’ respectively. The upper limit on the OFF:ON flux density is 
also denoted by Soff : Son - Note that we assume an upper limit of SNR = 3 for the OFF emission. The uncertainties in the flux densities are also quoted in 
parentheses with respect to the least significant digit. 


PSR 

V 

(MHz) 

Au 

(MHz) 

Non 

Noff 

Ton 

(s) 

Toff 

(s) 

SNR 

VFeq 

(ms) 

Son 

(mJy) 

Soff 

(/iJy) 

‘5'off • *5^011 
(10-3) 

Jf 634-5107 

f374 

288 

26 

136 

f0986.8 

85650.4 

I6f.O 

14.4 

0.40 (8) 

< 2.6(5) 

<7 (2) 


1518 

576 

13 

27 

7488.0 

f7072.6 

ff9.4 

15.7 

0.22 (4) 

< 3.6 (8) 

<17(5) 

JI7I7-4054 

732 

64 

2 

2 

874.9 

3784.9 

132.5 

69.6 

7(f) 

< 80 (20) 

< 11 (3) 


f374 

288 

44 

55 

f2I52.8 

16104.4 

f393.8 

13.3 

2.3 (5) 

< 4.4 (9) 

< 1.9(6) 


f5I8 

576 

8 

7 

2156.5 

2096.6 

406.5 

fl.6 

0.9 (2) 

<7(f) 

<7 (2) 


3094 

1024 

2 

2 

874.9 

3963.3 

158.2 

7.8 

0 .6(1) 

<5(f) 

<9 (3) 

Jf853+0505 

f374 

288 

27 

If 

16466.9 

3768.0 

305.8 

f05.7 

1.3 (3) 

< 27 (5) 

< 20 (6) 


f5I8 

576 

fO 

f 

8985.6 

898.6 

f90.9 

93.2 

0 .6(f) 

< 30 (6) 

< 50 (20) 


average flux densities in either the 10- or 50-cm band for these 
sources. 

Comparing the flux densities of PSRs J1634—5107, 
J1717—4054 and J1853+0505, we place upper limits on their 
Soft : Son values of < (7 ± 2) x 10"®, < (1.9 ± 0.6) x 10"® and 
< (2 ± 0.6) X 10“®, respectively. The stringent null confusion 
limit placed on PSR J1717—4054 indicates that the source could 
undergo deep nulls which are well below the sensitivity threshold 
of the Parkes observing system. Whereas, the limits placed on 
PSR J1634—5107 and particularly PSR J1853+0505 indicates that 
null confusion could occur in these sources. For PSR J1853+0505, 
this is consistent with its observed weak emission, which is 
detected just above the noise level (see also S I6.2t . 


4 TIMING 

Considering the F variations obs erved in several moding and/or 
nulling objects jLvne et al.ll2010h . we sought to characterise and 
compare the timing properties of the pulsars in our study. Here, 
we computed the best-fit timing solutions for the objects using the 
TEMP02 packag^ and the PSRCHIVE software suite. 

Due to the low number of 10/50cm ON observations, we 
only considered observations obtained at L-band for this analysis. 
Here, we formed analytic templates from the highest-SNR PAFB 
(all sources), PDFS (PSRs J1634-5107 and J1717-4054), LDFB 
(PSR J1853+0505) and LAFB (PSR J185 3+0505) profiles us¬ 
ing PA AS. Pulse time-of-arrivals (TOAs; e.g. [Manchester & TavloJ 
Il977h were generated through cross-correlation between the tem¬ 
plates and corresponding observations for each source using pat. 
Note that instrumental delays between receivers and/or backends 
were also fitted in tempo2. This process also ensures that delays 
caused by inaccurate alignment of the different templates used for 
the same source were accounted for. The results of this analysis are 
presented in Table[^and Fig.[3 

We find that the timing solutions for PSRs J163 4—5107 and 
J1853+0505 are consisten t with previous findings jHobbs et al.l 
1 20041 : UTorimer et ^l2006ll . In contrast to PSR J1853+0505, we 
note that the residuals of PSR J1634—5107 are not white, indicat¬ 
ing that variations may exist in this source. However, our observ¬ 
ing cadence is unsufficient here to infer anything conclusive about 
spin-down variations in the object. We also find substantial timing 


® An overview of this timing package is provided bv iHobbs et alj l2006h . 
See also http://www.atnf.csiro.au/research/pulsar/tempo2/ for more details. 
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Figure 5. Top panel: best-fit timing residuals for PSR J1634—5107, from 
an 11.1-yr span of observations, with crosses and squai'es denoting TOAs 
obtained at ^^1374 MHz and 1518 MHz, using the Multibeam and H- 
OH receivers, respectively. Centre panel: best-fit timing residuals for 
PSR J1717—4054, from 9.6 yr of data, with crosses and squares denot¬ 
ing TOAs obtained at 1374 MHz and 1518 MHz, using the Multibeam 
and H-OH receivers, respectively. Bottom panel: best-fit timing residu¬ 
als for PSR 11853+0505, from an 11.9-yr long data set, with circles, 
crosses and squares denoting TOAs obtained with the PAFB (1374 MHz 
and 1518 MHz), LAFB (1396 MHz) and LDFB (1402 MHz and 1520 MHz) 
backends, respectively. 
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Table 3. Fitted and derived properties of PSRs J 1634—5107, J1717—4054 and J1853H-0505 from their best-fit timing solutions. The distance of each source 
is taken from the NE2001 model, and is assumed to have a 20 % uncertainty. Errors are displayed in the parentheses, for parameters which were fitted, and 
are in units of the least-significant digit. 


PSR 

RA (J2000) 

(h: m: s) 

Dec. (J2000) 

u 

(s-1) 

u 

(10-1® s-2) 

DM 

(pc cni~^) 

d 

(kpc) 

Epoch 

(MJD) 

Vtoa 

Ts 

(yr) 

RMS 

ifis) 

J1634-5107 

16:34:04.99(8) 

-51:07:45.6(9) 

1.97100259922(5) 

-6.1167(5) 

373(2)“ 

6.1 

54420 

50 

11.1 

4537 

J1717-4054 

17:17:51.8(1) 

-41:03:20(5) 

1.12648304468(3) 

-4.6737(7) 

306.9(1)® 

4.7 

55035 

108 

9.6 

22011 

J1853-I-0505 

18:53:04.32(4) 

-1-05:05:29(1) 

1.10480469655(3) 

-1.5631(2) 

279(3r 

6.6 

53982 

89 

11.9 

7857 


“Lorimer et al. (2006), ^Kerr et al. (2014), ‘^Hobbs et al. (2004). 



Figure 6 . Variation in v against time for PSR J1717—4054, using fit win¬ 
dows of 200 d and stride steps of 50 d. The measured iy) from the overall 
timing solution is also overlaid. 


noise in PSR J1717-4054 (RMS - 22 ms: c.f. lHobbs et alj2odi) . 
which is comparab le to that observed in a large number of sources 
jHobbs et alEoiOh . 

The timing noise in PSR J1717—4054 can be interpreted as 
systematic v variations. To quantify this, we performed stride fits 
across i ts timing data usi ng tempo2, following the method out¬ 
lined in iLvne et al.l ( l2010h . Here, we used data windows of 200 d 
in length, and stride steps of 50 d, to provide the necessary timing 
accuracy while not compromising too heavily on time resolution. 
From this analysis, we find that the source exhibits a peak-to-peak 
spin-down variation Ai/pk = 1.53 ± 0.03 x 10“^® s~^ and a frac¬ 
tional spin-down variation | Ar'pk|/(f') = 33± 1 % over the course 
of hundreds of days (see Fig. [Hi. _ 

Following the method outlined in lYoung et alj(l2012h . we per¬ 
formed WWZ analysis on the v data to determine if the variations 
were periodic. However, we find that the cadence of our data is too 
poor to allow any periodic trend to be accurately determined. We 
were also unable to detect any significant pulse shape variation over 
time due to the narrow pulse longitude range over which emission 
is detected. This precludes any direct correlation between v varia¬ 
tion and pulse shape variation over time. 


5 POLARIMETRIC PROPERTIES 

We also analysed polarimetric Multibeam observations of 
PSRs J1634—5107 and J1717—4054. As no single pulse observa¬ 
tions were available in these data, we present only the average pulse 
polarisation properties. The pola risation calibration was carried out 
following the scheme outlined in IWeltevrede & Johnston! l l2008h . 


5.1 PSR J1634-5107 


We used 12 Multibeam observations to probe the polarimetric prop¬ 
erties of PSR J1634—5107. These were aligned, using the timing 
solution presented in Table and then averaged to produce an 
integrated profile from the 42 min of data using the PSRCHIVE 
software suite. We determ ined the rotation measure (RM; see, e.g., 
iLorimer & Krameil2005ll of the source. RM = —151±8 rad m~^, 
for the first time using the rmfit package jNoutsos et alfcOOsh . 

After correcting for RM, we analysed the intrinsic Stokes pa¬ 
rameters (I,Q,U,V) and degree of linear {L/I = /1) 

and circular (V/I) polarisation for the source (see top panel of 
Fig.jT}. Using a 3cr limit for the polarisation data, we found that the 
pulsar exhibits high linear ((L/7) = 19 ± 2 %) and modest circular 
((|U|//) = 10±1 %) polarisation, when averaged over bins within 
the 10 % intensity pulse width (VFio) range. We also computed the 
maxima of the polarisation profiles; i.e. (I///)max = 40 ± 10% 
and (|U|//)max = 30 ± 20 %. Here, we quote uncertainties on the 
degrees of linear and circular polarisation in terms of their respec¬ 
tive quadrature errors; e.g. AL = where ctl is the RMS 

of L in the off-pulse region and Non is the number of on-pulse bins 
considered). 

Further to the above, we sought to determine the magnetic in¬ 
clination angle of the source (a) and the impact parameter of the 
LOS (/3). For this analysis, we fitted the polarisation position an¬ 
gles (P As) of the integrated emission wi th the rotating vector model 
(RVM; iRadhakrishnan & Cookd 1969I). adopt i ng the minimi¬ 
sation technique presented by Rookvard et al.l ( 1201 5ll to obtain fit 
constraints with 3a PA values. 

We find that the best-fit to our data is consistent with the RVM 
(see Fig. |7j. However, we cannot place rigorous constraints on 
the a and P parameters due to the low number of significant PA 
data points available. Instead, if we assume a = 90°, we can es¬ 
timate the maximum P valu e (/Sqq) from the maximum gradient 
(dT'/ci(/))max of the RVM fit ( lKomesarofllll970l : see below). In the 
general case, P is less than this value: 


Pgo = sin 


sin (90°) ( ^ 


dp 


( 1 ) 


5.2 PSR J17I7-4054 

We used 35 Multibeam observations to characterise the polarimet¬ 
ric properties of PSR J1717—4054. The 146 min of data were 
aligned and folded, using the same method as in § [511 to produce 
an integrated profile for this analysis. We again used rmfit to mea¬ 
sure an RM = —809 ± 3 rad m~^ for the source. This best-fit 
RM va lue, which represents a refinement of the result of iKerr et al.l 
( |2014|) . was then used to correct the average profile, resulting in the 
profile shown in the top panel of Fig.j^ 
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Figure 7. The best RVM fit to the integrated 42 min observation of 
PSR J1634—5107. Top panel: the average emission properties of the inte¬ 
grated observation, showing the total intensity profile, as well as the lin¬ 
ear and circular polarisation profiles. Bottom panel: integrated PA data 
(crosses), with the best-fitting RVM (solid line) and location of the PA 
swing inflection point (ijio; dotted line) overlaid. 



0 50 100 150 200 250 300 350 


Pulse longitude (deg.) 

Figure 8 . The average emission properties of the 146 min observation of 
PSR J1717—4054. Top panel: the total intensity (black line), linear (red 
line) and circular (blue line) polaiisation profiles. Bottom panel: integrated 
PA data for the average profile. 

We find that the source exhibits modest linear ((L/7) = 8.1 ± 
0.1 %) and circular polarisation {{\V\/1) = 4.9 ± 0.1 %), when 
averaged over the Wio range (with (Z//I)max = 30 ± 4% and 
(14//)max = 6.9 ± 0.2 %). We also attempted to fit the PA data for 
the source, using the same method as in S I5.1I but could not obtain 
any constraints on the values of a or from the resultant reduced 
surface. This is due to the clustering of the PA data points which 
are not highly constrained by the RVM (see bottom panel Fig.[^. 


6 DISCUSSION AND CONCLUSIONS 
6.1 Emission Variability 

Thorough analysis of our extensive data set has shown that 
PSRs J1634-5107, J1717-4054 and J1853+0505 exhibit long 
term nulls and emission phases (i.e. minutes to many hours), 
as well as substantial NFs (i.e. ~ 67 % — 90 %). Remarkably, 
PSR J1853+0505 was also shown to exhibit a weak emission state, 
in addition to its strong and null states (see S I6.2t . 

Short timescale nulling (~ 1 — 40 P) was discovered in 


PSR J1717—4054, during its ON phases, which acts to modulate 
the NF and average detection rate of the source over sub-integration 
timescales. This behav iour was al s o disc overed in an independent 
analysis performed by iKerr et aP (|2014|) . who analysed the same 
archival data and two dedicated 9.5-h single-pulse observations. 
The latter observing runs constrained a number of mode transi¬ 
tions and thus allowed tentative limits to be placed on the aver¬ 
age ON and OFF emission timescales of ton ~ 1200 ± 700 s and 
ton = 7000 ± 5000 s, respectively. These average emission du¬ 
rations are found to be consistent with the results obtained in this 
work, and indicate that the source exhibits a wide range of emission 
timescales. Thus accounting for varied NF reports in the literature. 

Overall, the above results place PSRs J1634—5107, 
J1717—4054 and J1853+0505 at the extreme end of the ‘nulling 
continuu m’, where there is a current deficit of know n objects 
(see, e.g., iKeane et ^1201 ll : ISurke-Spolaor et aOlTOl ih . For such 
sources, it is clear that long observations (> 1 h) are required 
to best character ise their emission properties, as demonstrated in 
lKerretal.U2014h . 


6.2 Deep Nulls or Weaker Emission States? 

PSR J1853+0505 is shown to exhibit three emission states: a weak, 
strong and null state. During the weak mode, emission is barely de¬ 
tected above the noise level in our observations. This behaviour is 
remarkably similar to that observed in PSR J1107—5907, where 
emission at the lowest end of its w eak-mode pulse-en ergy distri¬ 
bution is easily confused with nulls jYoung et al.ll2014l) . Similar to 
PSR J1107—5907, integration of > 10^ pulses is required to detect 
the particularly weak emission of PSR J1853+0505, leading to an 
average SoB ■ -Sweak ^ 0.11 ± 0.03 at 1374 MHz. These simi¬ 
larities indicate a close connection between the two objects which, 
in turn, advocates further, more detailed analysis of the different 
emission modes of PSR J1853+0505. 

For the other sources in our study, PSRs J1634—5107 and 
J1717—4054, only two discrete modes were observed in each ob¬ 
ject. From this work, we place null confusion limits of < (7 ± 2) x 
10“^ and < (1.9 + 0.6) x 10“^ on these pulsars, respectively. Con¬ 
sidering the fact that PSR J1634—5107 is approximately 3 times 
fainter than PSR J1853+0505, it is unsurprising that observations 
of the former did not result in the discovery of a weaker emission 
state. As such, PSR J1634—5107 could exhibit a particularly weak 
emission state which is not probed by our observations, or it could 
actually undergo deep nulls. In the case of PSR J1717—4054, we 
suspect that the source exhibits t he latter. This is supported by in¬ 
dependent analysis performed bv iKerr et al.l (|2014|) . which resulted 
in a more stringent limit on its null confusion limit (i.e. < 10“^). 
Thus, our observations clearly do not probe a potential lower inten¬ 
sity emission state in this object. 

Since particularly weak emission has been reported in this 
work, and in a number of o ther pulsars (e.g. 5weak : -Sstrong 
2 X 10“^ for J1107— 1107: lYoung et al.ll2014l) . it is possible that 
pulse nulling may only represent an instrumental sensitivity bias 
on certain members of the pulsar population. Overall, these findings 
provide additional motivation to perform more regular and longer 
observations of nulling objects, which will likely only be possible 
with next generation telescopes such as MeerKAT and the SKA. 
Such steps would ultimately help to construct a census of nulls and 
mode transitions in the pulsar population, which is required to help 
form a fully consistent pulsar emission model. 
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6.3 Timing Properties 


REFERENCES 


We have confirmed and updated the timing solutions of 
PSRs J1634-5107, J1717-4054 and J1853+0505, using the 
longer data spans available in this work. Both PSRs J1634—5107 
and J1853+0505 are found to exhibit weak timing noise, in spite 
of their extreme variability. Whereas, PSR J1717—4054 is shown to 
exhibit modest timing noise, similar to that observed in a number of 
pulsars which switch b etween magnetospheric states dptobbs et al.l 
l20ld:lLvne erZll2010h . 

Upon further investigation of the timing behaviour of 
PSR J1717—4054, we found that its timing noise can be ac- 
counted for by su bstantial variation in v over hundreds of days (c.f. 
lUvne et aDl2010h . This variation is too slow to be associated with 
transitions between ON and OFF phases in its emission. Therefore, 
it is likely that the source also exhibits emission variability over 
similarly long timescales. While we were unable to confirm this 
hypothesis through analysis of our data, we predict that more fre¬ 
quent detections of this source with high-time resolution data will 
lead to a direct association with z> variation. Overall, it is clear that 
PSR J1717—4054 is a multi-state switching object, which can be 
used to probe both short- and long-timescale variations in pulsar 
magnetospheres. 


6.4 Polarimetric Findings 


We have presented the first measurement of the RM for 
PSR J1634—5107, which has enabled us to fit its a and ft parame¬ 
ters. While we were unable to place constraints on a, we were able 
to infer Pqo ~ 3°. Under the assumption that radio beams are com¬ 
prised of core and conal emission components, we can determine a 
from the observed pulse width of the core fFFfn'^‘^: lRan ^Il990ll : 


a = sin 


/ 2.45° \ 

VFF5°o°"'°PO V ■ 


( 2 ) 


The emission profile of J1634—5107, can be described by a strong 
core component, flanked by two conal emission components (see 
Fig.|7j. We thus estimate = 6.9° and a ~ 30°. Flowever, 

we stress that longer polarimetric observations of the source in its 
ON state, with higher time resolution, are required to better con¬ 
strain its emission geometry and verify our results. 

We also estimate an RM = —809 ± 3 rad m“^ for 
PSR J1717—405 4. This value is consistent with the result of 
I Kerr et aP ( 1201 4h and is found to be one of the highest values 
measured for a puls ar not located in a globular cluster (see, e.g., 
iNoutsos et al.ll20()8 and references therein). After correcting the 
source’s integrated profile for RM, we also tried to fit the RVM 
to its PA data. Flowever, we were unable to obtain any constraints 
on a or P due to the clustering of PA data points across the narrow 

S iulse p rofile, which is again consistent with the result of iKerr et al.l 
20141) . Nevertheless, further investigation into the polarimetric 
properties of the source would be interesting to shed light on its 
short-timescale variability. 
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